Cooper pair cotunneling in single charge transistors with dissipative electromagnetic 

environment 
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We observed current-voltage characteristics of superconducting single charge transistors with on- 
chip resistors of R ~ Rq = h/4e 2 ~ 6.45 kf2, which are explained in terms of Cooper-pair cotunnel- 
ing. Both the effective strength of Josephson coupling and the cotunneling current are modulated 
by the gate-induced charge on the transistor island. For increasing values of the resistance R we 
found the Cooper pair current at small transport voltages to be dramatically suppressed. 

PACS numbers: 74.50. +r, 74.25.Fy, 74.40. +k 



During the last decade superconducting circuits with 
mesoscopic tunnel junctions have attracted much at- 
tention of researchers. Remarkable quantum properties 
of these circuits, which stem from interaction between 
Josephson and charging effects, enable the construction 
of charge-based qubits , sensitive electrometers Q , am- 
plifiers |3j], etc. Since the Josephson supercurrent itself 
is of an elastic nature, tunneling of Cooper pairs in such 
circuits is dramatically influenced by the Coulomb inter- 
action and dissipation. 

The tunneling in the two-junction system (single 
Cooper pair transistor) with a dissipative environment 
was recently theoretically studied by Wilhelm et al. 3 
and experimentally investigated by Kycia et al. [j| and 
Lu et al. @ who used the advantage of locally tunable 
dissipation associated with the conductance of 2DEG be- 
neath the transistor structure. These papers were fo- 
cused on the regime of sequential tunneling of Cooper 
pairs (STCP) dominating when the electric potential of 
the transistor island was tuned "in resonance" (see pa- 
rameter diagram in Fig. 3 of Ref. Q). Although a quan- 
tum regime of simultaneous pair tunneling across both 
junctions, i.e. cotunneling of Cooper pairs (CTCP), was 
roughly evaluated in Ref. Q, it had not been, to our 
knowledge, studied experimentally. 

Under the term "cotunneling" we understand discrete 
tunneling of single Cooper pairs across a whole tran- 
sistor [7J. Although numerous experiments show that 
Cooper pair transistors without deliberately attached re- 
sistors exhibit either a gate-dependent switching current 
at ramp-current bias or a current peak at low voltage 
bias (see, e.g., Ref.Q), these data indicate to a through 
supercurrent. Due to unavoidable dissipation leading to 
fluctuations of the overall Josephson phase <j>, the ob- 
servable value of the supercurrent is non-zero at finite 
transport voltage. Phase <j> is still a good variable, while 
the number of transferred Cooper pairs is meaningless, 
so the term "single" is not applicable to this regime. 

When an appreciable impedance is introduced into the 
network, i.e. ReZ = R ~ Rq = h/Ae 2 « 6.45 Ml, 
quantum fluctuations of phase <f> increase dramatically, 
and this leads to incoherent tunneling of individual pairs 



across the transistor, i.e. CTCP. In this paper we investi- 
gate the CTCP effect in Al transistors with thin-film Cr 
local resistors of different values R inserted into the bias 
leads. We observe the gate and bias voltage dependencies 
of CTCP current and demonstrate that a significant sup- 
pression of CTCP at small bias voltage can be achieved 
by a further increase in the resistance R. 

The Cooper pair transistors were fabricated using the 
angle evaporation technique for three metal shadows 
(Cr/Al-1 (oxidized)/ Al-2) evaporated at different angles 
through the same mask. Details of the fabrication pro- 
cess can be found elsewhere . We put on the same chip 
a number of nominally identical symmetrical transistors 
(with tunnel junctions of type A1/A10 X /A1 with dimen- 
sions 40 nm x 80 nm) differing only in the lengths of the 
Cr microstrips (each 0.1 /jm - wide and 7 nm - thick), 
ranging from 0.32 fj,m to 3 //m. The specific film resistiv- 
ity was from 5.5-7 Ml per /zm of length yielding values 
of resistance R in the range of 2-20 Ml. 

The measurements were performed in the dilution re- 
frigerator at the base temperature of about 20 mK. The 
sample was mounted in a copper shielding box. All bias 
and gate lines were supplied with microwave filters, made 
of Thermocoax® cable located in the coldest part of the 
setup. The electric diagram of the sample bias is shown 
in the upper inset on Fig. 1. The normal-state tun- 
nel resistance of the transistors, Rn = Ri + R%, where 
i?i.2 are tunnel resistances of individual junctions, was 
remarkably uniform and ranged from 40 Ml to 44 Ml. 
The superconducting gap and the corresponding values 
of the Ambegaokar-Baratoff critical current of the junc- 
tions (assuming symmetry of the transistors) were es- 
timated to be Aai ~ 200 /xeV and 1° ~ 16 nA, respec- 
tively. These data yield a characteristic Josephson energy 

30 /ieV, where 



of junctions E J1:J2 = Ej = §£/° 



$o = 2 * s the flux quantum. From the normal-state 
I— V curves we evaluated a characteristic charging energy 
of the transistor island E c = e 2 /2Cs rs 150-160 /ieV, 
where Cs = G\ + G<x + C g ~ 530 aF (typical value) is the 
total capacitance of the island and C1.2 = C « 260 aF 
and C g ss 9 aF <C C are the capacitances of junctions 
and gate, respectively. Due to the small C g , the gate- 
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induced charge is fixed, Qq — C g V g + Qoo, where Qoo 
is the background charge. The reasonably small ratio of 
characteristic energies in our samples, A = Ej/E c as 0.2, 
allowed us both clear observation and transparent inter- 
pretation of the CTCP effect (<x A<J). Mor eover, linearity 
and small resistor lengths let us consider them as lumped 
elements (pure resistors) and neglect in our model their 
distributed self-capacitance (ss 60 aF/^m). 

I — V characteristics of superconducting single charge 
transistors are generally featured by several contribu- 
tions to the net current. These contributions result 
from the tunneling of Cooper pairs, single quasiparti- 
cles and combining tunneling of both pairs and quasi- 
particles. At small transport voltage and temperature 
(eV^ksT >C 2_E c ,2Aai), the tunneling of quasiparticles 
cannot make a notable contribution to the net current 
because of the exponentially large subgap resistance. In 
this limit, the Cooper pair (Josephson) current compo- 
nent is dominant. Still, the rare events of quasiparticle 
tunneling (poisoning of the island) at \Qq\ > e/2 were 
possible and they resulted in le-periodic gate dependence 
of the I — V curves (see bottom inset in Fig. 1). There- 
fore, slow ramping of the gate voltage (when measuring 
gate-modulation curves) provided cyclic variations of the 
polarization charge on the island Qq in the range from 
—e/2 to e/2. Within this range of Qq, change of the 
island charge by 2e would require an energy SE c h rang- 
ing between 2E C (at Qq = ±e/2) and 4E C (at Qq = 0). 
Therefore, the process of Cooper pair transfer with "land- 
ing" on the island, i.e. the STCP, was exponentially 
suppressed, tx exp(— SE^/hsT), so only the CTCP con- 
tributed to the net current. 

Figure 1 shows the I — V curves of three samples 
with resistance R of nominally zero value (the actual 
impedance in this case is of the order of the vacuum 
impedance Zq rj 377 fl, although its frequency depen- 
dence is unknown), 3 kS7 and 5 kfi. One can see a re- 
markable change from the peak (cf. Fig. 1 in Ref. 0) to 
step-like and then to an almost linear shape of the I —V 
curve that can be interpreted as a growth of fluctuations 
of the total Josephson phase leading to a depression of 
the supercurrent feature. 

A similar behavior of a single junction was described 
by Averin et al. |10| in terms of incoherent tunneling of 
pairs yielding current I — 2e[T(2eV) — T(—2eV)] and a 
tunneling rate 

T{E) = ^h E2jP{E) > (1) 

where the function P(E) with E = 2eV describes the 
probability of energy dissipation in the environment net- 
work. The experimental observation of this regime was 
reported by Kuzmin et al. [ill ]. 

In the case of a transistor, one may consider Ej in 
Eq. (JTJ as an effective Josephson coupling dependent on 
the gate voltage V g . To support this idea, we observed 
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FIG. 1: The electric circuit diagram (upper inset) and I — V 
curves of three Cooper pair transistors with different nomi- 
nal values of the resistance R (circles - 0, squares - 3 kfi and 
triangles - 5 kfi) at gate voltage values V g corresponding to 
maximum of current. The solid line shows computed data 
assuming a symmetric transistor with E.j = 26 /ieV and a re- 
sistor of R = 3.3 kf2 with an effective temperature of 90 mK. 
The gate dependence in the transistor with R — 3 kf2 at 
V — 30 jitV is presented in the bottom inset, where the solid 
(dashed) line shows the data calculated within an experimen- 
tally accessible (inaccessible) range. 



the shapes of the I — V curves similar to those given by 
the P-theory for single junctions with ohmic environment 
[ToL Il2| . Moreover, the observed shapes were immutable 
to V g whose variation led only to scaling of the current. 

To model the behavior of our I —V curves, we present 
the Josephson coupling term of the transistor Hamilto- 
nian, Hj = —Ej± cost^i — Ej 2 cos</?2, hi the form 

Hjfatp) = -£#(0cos[¥> + 7(0]. (2) 

where Ej = (E 2 n + E% + 2E J1 E J2 cos (f) 1/2 , tan 7 = 
[{Eji -E J2 )/ {Eji + E J2 )} tan(0/2) and 4> = (p x + ip 2 
(ipi is the phase drop across the z-th junction). The 
difference phase ip = \{<£i — (fi 2 ) is a variable conju- 
gate to the island charge Q — —2ei-j^. Q enters the 
charging energy term H ch (Q) = E C [(Q + Qo)/e} 2 . For 
a given value of operator (j), the transistor Hamiltonian 
H tI = H c h(Q) + Hj((f>,<p) is similar to that of a single 
junction. Taking the fact into account that the com- 
mutator \<j), i?tr] = 0, we can diagonalize Hti(<p, <p, Q) = 
i?tr(</>, q, n) in the basis of the Bloch states \q, n), where q 
is the quasicharge and n = 0, 1, 2, ... the band index |l3j |. 
Assuming that only the ground state, n = 0, is occupied 
(i.e., the effects of thermal, photon-assisted and Zener 
breakthrough excitations of higher bands are negligibly 
small) and noting that variable q is decoupled from envi- 
ronmental degrees of freedom due to small C g and fixed 
by V g , we present the transistor Hamiltonian as 

H %v = E [\{^),Q ]. (3) 
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Here Eq is the ground state energy of an equivalent sin- 
gle junction with variable ratio A = Ey((j))/E c < 2\q 
and q — Qo- E is 27r-periodic function of <fi and 2e- 
periodic function of charge Qo and its shape depends on 
the relation between Eji t j2 and E c . 

Due to the periodic dependence on 0, the Hamiltonian 
Eq. can be expanded into a Fourier series 
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(4) 



fc=0 



For Qo 7^ e mod (2e) or/and at notable asymmetry of 
the transistor, Ej\ ^ Ej2, the terms with k = and 1 
are dominating in Eq. (@J . This is in accordance with al- 
most sinusoidal supercurrent-phase relation 0] with ef- 
fective critical current I' r (Qo) = gr-^j (Q ) modulated 
by gate (see corresponding plots for different values of A 
in Fig. 2a of Ref.|l5j). In particular, in the symmetric 
transistor with Ao -C 1, i' r can be found explicitly from 
Eqs. (4)- (5) of Ref. using the wave functions given 
by Eq. (All) of Ref. Ha, i.e., 



I t c I (Qo) = ^ I- ((),:<) 



1° 



(5) 



In the degenerate case of Qo — e mod (2e) and Ej\ — 
Ej2 = Ej < E c , the ground state energy Eq. © is E ~ 
— Ej | cos ^ | +const, while the maximum critical current 
is /* r (e) « 0.5/ c i iC2 . The Fourier coefficients in Eq.QJ 
are equal to 



Ef\e) = - 



4 (-1) 



fc-i 



7T {2kf 



-Ej 



k > 1, 



(0) 



so the amplitudes of the principal (first harmonic) term 



E { j\e) w 0A2Ej and the ratio \E ( f > (e)/E ( J 1> (e)\ = 0.2. 

The expansion Eq. (0} presents the transistor Hamil- 
tonian as a weighted sum of operators e ±lk ^, each de- 
scribing tunneling of a cluster of k Cooper pairs through 
the transistor in the positive (negative) direction. This 
presentation enables particular rates for tunneling of k 
pairs, T k (E), to be found and the total current to be 
calculated via the corresponding sum. 

Due to finite impedance R, the phase <p couples to 
the environmental degrees of freedom via an operator of 
phase on the resistor <^r: cf) — ^Vt — <pR. Phase 4>r is 
related to the infinite set of bath variables. The effect 
of this coupling can be treated by means of function P 
determined by the charge of the tunneling particle (clus- 
ter of particles), temperature and the real part of total 
impedance Z as seen by the tunneling charge fl6| . 

In our case, the real part of the impedance is equal to 



(2), 



ReZ(w) = R 1 + (ujRC tr ) 



(7) 



where the transistor capacitance Cor- 
responding to the characteristic charging energy of 4_E C ) 



C1+C2 2 ^ LU1 



< 

Q. 



Fit: 1 = Ax V 2p ~ 1 

P = 0.86 A^fe^ 
R fit = 5.5k£l j^sr 




^ p = 1 .5 c<y* 
R fit = 9.7kn / p = 

■ / V 


2.62 J 
>.9 knSv 
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FIG. 2: The I — V curves of three Cooper pair transistors with 
different values of resistance R (up-triangles - 5 kf2, diamonds 
- 9 kfi and down-triangles - 17 kf2, cf. the fit values Rat) and 
their power-function fit in accordance with Eq. 1101 . The 
fitted values of prefactor A = IoVq~ 2p were used to restore 
the values of Ej (= 33, 34 and 28 peV, respectively). 



and the tunneling charge is equal to 2fce. The rates, 
T k {E) = ^-[Ef ] ] 2 P(E), and the total cotunneling cur- 
rent, 



I cot = 2eY,k [T k {2keV) - T k {-2keV)}, 



k=l 



are found perturbatively, assuming 0.5Ej / 'AE C = Ao/8 <C 
1. If Qo not close to e, the term with fc = 1 is essential, 
and the cotunneling current is given by the formula sim- 
ilar to that derived by Ingold et al. [12j for tunneling 
across a single junction at 2eV <C h/RC, 

loot- e p 



2E C 



\T[p-j(eV/TTk B T)}\' . / eV 
x : sinn 



r(2 P ) 



knT 



, (9) 



where p — R/Rq, T is Gamma function and 7 rs 
0.577 Euler's constant. The perturbation theory result, 
Eqs. ©-©, is valid if 7 cot < p/* r . The latter condition 
means that the cotunneling events are rare and incoher- 
ent and this is nicely obeyed in our experiment. 

Expression Q pretty well gives the transport and 
gate voltage dependencies as well as the magnitude of 
the measured current (see calculated curves in Fig. 1). 
The values of the fitting parameters R and Ej are in 
good agreement with their estimated values (10% and 
13% deviation, respectively), while the effective temper- 
ature T e g of the resistor was found to be 90 mK (cf. 
bath temperature of 20 mK). This value matches with 
T c ff = 70 — 80 mK found in the same setup in measure- 
ments of the three- jun ction Al-Cr single electron pump 
with local resistors [17J and can be ascribed mostly to a 
noise effect. 

We achieved the CTCP regime in the range of the gate 
charge \Qo\ < f ■ The gate modulation depth of cotun- 
neling current was found to be d = / C ot(§)/7 C ot(0) ~ 2, 
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that is in agreement with the value given by Eq. JSJ, 
i.e. d = [^ r (§)//* r (0)] 2 = f ~ 1.78 0. Approaching 
the points Q = ±e one could expect an increase in the 
cotunneling current by a factor of (2/Ao) 2 ~ 10 2 , but 
in these points STCP dominates because effective dis- 
sipation for tunneling across either junction is notably 
smaller than that for CTCP, p ^ p/4 |. As a mani- 
festation of such behavior in our samples, we observed a 
current increase by more than two orders of magnitude 
when we applied a fast (~ 1 MHz) gate voltage sweep of 
amplitude > e/C g . Such sweep allowed quasiparticle poi- 
soning of the island to be avoided for short time intervals 
when passing the regions Q « ±e. 

Finally, an issue of particular interest is the suppres- 
sion of CTCP at larger dissipation, i.e. at R > Rq 
The increase in R should lead to a blockade behavior 
of the I — V curves at small V and development of a 
broad current peak moving to V = The latter fea- 
ture can hardly be clear observed in experiment because 
at such a large voltage other current components (quasi- 
particle and combined pair-quasiparticle tunneling) are 
significant ^lj. On the other hand, for small voltages 
between ksT/2e and e/C, cotunneling still dominates 
and for its estimation one can apply the zero tempera- 
ture expansion of function P and will get 



Io 



V~o) 



where In 



7rpe- 2 T$ (^ r ) ; 
32T(2p)E c 



(10) 



and Vq = Thus, the CTCP current obeys the power 
law, i.e. it is efficiently suppressed at large p. This ef- 
fect was demonstrated and the results are presented in 
Fig. 2. Although the experimental voltage range was 
rather limited, the data are in good agreement with the 
dependence Eq. (fTUf) for parameter values close to those 
estimated from measurements. 

In conclusion, the obtained results on measurements of 
CTCP are also of practical importance, e.g., for improv- 
ing the performance characteristics of the Cooper pair 
pumps. In the latter devices, CTCP may significantly 
degrade the accuracy of charge pumping [2(| ■ In contrast 
to the normal junction case, an increase in the number 
N of junctions in the superconducting pump array seems 
to be not suitable for suppressing unwanted CTCP [2l| . 
On the other hand, the use of local resistors 0, H3| en- 
sures such suppression and, therefore, is promising for 
the construction of an accurate Cooper pair pump with 
a reasonable number (N = 3 — 5) of junctions in the 
array. 

This work was partially supported by the EU through 
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